Disordered quantum walk-induced localization of a Bose-Einstein condensate 
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We present an approach to induce localization of a Bose-Einstein condensate in a one-dimensional 
lattice under the influence of unitary quantum walk evolution using disordered quantum coin opera- 
tion. We introduce a discrete-time quantum walk model in which the interference effect is modified 
to diffuse or strongly localize the probability distribution of the particle by assigning a different set 
of coin parameters picked randomly for each step of the walk, respectively. Spatial localization of the 
particle/state is explained by comparing the variance of the probability distribution of the quantum 
walk in position space using disordered coin operation to that of the walk using an identical coin 
operation for each step. Due to the high degree of control over quantum coin operation and most 
of the system parameters, ultracold atoms in an optical lattice offer opportunities to implement a 
disordered quantum walk that is unitary and induces localization. Here we present a scheme to use a 
Bose-Einstein condensate that can be evolved to the superposition of its internal states in an optical 
lattice and control the dynamics of atoms to observe localization. This approach can be adopted to 
any other physical system in which controlled disordered quantum walk can be implemented. 



I. INTRODUCTION 

Localization of waves in disordered media originally 
was predicted in the context of transport of electrons 
in disordered crystals by Anderson Anderson local- 
ization, resulting in the absence of diffusion, originates 
from the interference between multiple scattering paths 
(cf. 0). This phenomenon is now ubiquitous in physics 
[H, HI ; it has been experimentally observed and theoreti- 
cally studied in a variety of sy stems, including light waves 
|5Hlfj| and matter waves fllHl4|. 

Quantum walks (QWs) (l5Hl9| . which are the quan- 
tum analog of the classical random walks (CRWs), evolve 
in position space involving interference of amplitudes of 
multiple traversing paths. The quantum features of inter- 
ference and superposition are know to make probability 
of QW spread quadratically faster with time than its clas- 
sical counterpart in one- dimension. Some studies have 
shown the localization of the interference of amplitudes 
between multiple traversing paths of the QW distribution 
around the origin from various different perspectives [20- 
26| . In particular, it is shown that the localization of the 
walk dynamics in one dimension can be controlled by in- 
troducing drifts with constant momentum between two 
consecutive steps of QW [20[ or by evolving the walk in a 
random medium characterized by a static disorder (22[. 

The key factor for the interference effect to result in 
localization is the broken periodicity in the dynamics of 
the system induced by the disordered media. However, 
broken periodicity need not be mediated by a disordered 
or a random medium alone; operations dening the dy- 
namics of the system can be made random to break the 
periodicity such that they mimic the effect of a random 
medium in the system and manifest localization. Taking 
this into consideration, we can carefully construct a QW 
evolution on a physical system such that the dynamics 
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of the walk without any disorder in the lattice is simi- 
lar to the dynamics of a walk in a disordered medium 
or in the lattices leading to localization. Owing to the 
high degree of control over most of the system parame- 
ters and recent experimental implementation of QW [13] , 
ultracold atoms in an optical lattice offer opportunities 
for the study of disordered QW-induced localization. 

Using the exibility and control that has been achieved 
over the ultracold atoms, Bose-Einstein Condensate 
(BEC) has been used for the study of disorder- induced 
localization. Using a cigar-shaped noninteracting BEC, 
exponential localization [12j and a crossover between ex- 
tended and exponentially localized states [Rlhave been 
experimentally demonstrated. Billy et al. [l2| demon- 
strated exponential localization of 87 Rb atoms when re- 
leased into a one-dimensional waveguide in the presence 
of a controlled disorder created by a laser speckle. Roati 
et al. [13], using 39 K atoms, demonstrated its localization 
in a one-dimensional bichromatic optical-lattice potential 
created by the superposition of two standing- wave polar- 
ized laser beams with different wavelengths. In the pre- 
ceding two experiments, localization was demonstrated 
through investigations of the transport properties and 
spatial and momentum distributions. Numerical study 
of localization of BEC in a bichromatic optical-lattice 
potential [28j and in a random potential [29j has been 
reported. Recently, it was shown that the localization 
by bichromatic potentials is produced by a trapping by 
the potential and is not due to a quantum suppression, 
in contrast to the Anderson model (30[. There are also 
other time-dependent phenomena that can induce popu- 
lation imbalance of two self-interacting BEC and, hence, 
localization in one of the double-well potentials unlike 
Anderson localization, which is a stationary phenomenon 
31]. 

Taking the preceding points into consideration, in this 
paper, we present a new scheme to observe dynamic lo- 
calization of ultracold atoms in an optical lattice. Peri- 
odicity in the dynamics of atoms in an optical lattice is 
broken using a disordered evolution of the QW and this 
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leads to the interference of amplitude of multiple travers- 
ing paths of atoms in an optical lattice to localize. Direct 
control over the quantum coin operation makes it possible 
to choose the random set of coin operations and control 
the dynamics of the QW [l8, 32], which, in turn, allows 
us to break the periodicity and localize the evolution. 
In particular, we discuss localization using a BEC with 
noninteracting and interacting atoms, respectively. This 
scheme can be expanded to any of the physical systems 
on which the QW can be implemented. 

This article is arranged as follows. In Sec. [Til we de- 
scribe the DTQW model on a line. In Sec. IIII1 we in- 
troduce the disordered DTQW model using controllable 
quantum coin operation randomly picked for each step of 
walk and present both diffusive and localization effects, 
respectively. In Sec. IIV1 we discuss the implementation 
of QW in ultracold atoms and localization of BEC. Fi- 
nally, we conclude in Sec. [V] 



II. QUANTUM WALK 

Like their classical counterparts, QWs are also widely 
studied in two forms, continuous-time QW (CTQW) |19j | 
and discrete-time QW (DTQW) [H EE [M H3 , and are 
found to be very useful from the perspective of quan- 
tum algorithms (35l - [38j . Furthermore, they have been 
used to demonstrate the coherent quantum control over 
atoms [39| and to explain phenomena such as the break- 
down of an electric-field driven system [2l| and wavelike 
energy transfer within photosynthetic systems [10, l4lj |. 
Generation of entanglement between two spatially sepa- 
rated systems is another application of QWs (26| . Exper- 
imental implementation of QWs has been reported with 
samples in nuclear magnetic resonance (NMR) systems 
[42l [43l ] , in the continuous tunneling of light fields through 
waveguide lattices 44j , in the phase space of trapped ions 
[45l based on the scheme proposed by [47| , with sin- 
gle optically trapped atoms (27|, and with single photons 
[48l |49| . Recently, implementation of a QW-based search 
algorithm in a NMR system has been reported [50] . Var- 
ious other schemes have been proposed for their physical 
realization in different physical systems [5l|-[54j . 

DTQW, the dynamics f which can be controlled by 
controlling the quantum coin operation, is used for the 
study presented in this paper. It is modeled as a particle 
consisting of a two-level coin (a qubit) in the Hilbert 
space H c , spanned by |0) and |1), and a position degree 
of freedom in the Hilbert space Ti p , spanned by \ip x ), 
where x £ I, the set of integers. A t-step DTQW with 
unit time required for each step of walk is generated by 
iteratively applying a unitary operation W, which acts 
on the Hilbert space H c <8> "H p : 



|**> =W*\9in), 



(1) 



is an arbitrary initial state of the particle at the origin 
and W = S(B ® 1), where 

B - B ^ = { e"« sin(0) -e"* cos(0) ) U{2) 

(3) 

is the quantum coin operation. S is controlled-shift op- 
eration 



hk-iMM + |i>(i|® h^ + i>ftM 



( 4 ) 

The probability to find the particle at site x after t steps 
is given by p(x,t) = (ip x \tr c (\^ t )(^ t \)\i; x ) . 
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FIG. 1: (Color online) Probability distribution of the QW 
evolved by assigning quantum coin operation in a U(2) group 
with random parameters to each step of the walk. Pa- 
rameters are assigned from £,9,( £ {0, 7r/2} for B^^,( an d 
Bo t e,o respectively. The distribution is after 100 steps of the 
walk on a particle with initial state at the origin \^i„s) = 
^7f(|0) -H|l}) g> \ipo)- We see that the variance of the distribu- 
tion is very much close to the variance of the classical random 
walk (CRW) distribution. 



III. DISORDERED DISCRETE-TIME 
QUANTUM WALK 

Direct control over the quantum coin operation B 
makes it possible to control the dynamics of the DTQW 
[3 HH ■ For a walk evolution on a particle in one dimen- 
sion with initial state 
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where 



at the origin using an unbiased coin operation, that is, 
Be £ with £ = £ = 0, the variance after t steps of walk is 
(cos((5/2)|0) + sin((5/2)e i0 |l)) ® |i/> ) (2) [l-sin(6')]i 2 [33] and a symmetric probability distribution 
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in position space is obtained [3J, l55[ . Nonzero values for 
parameter £ and C when £ ^ C introduces asymmetry to 
the probability distribution [32j. It is also shown that 
DTQW with S ^ n/2 in Eq. @ returns asymmetric 
probability distribution, and parameters £ and C can be 
adjusted to make the distribution symmetric. 

In a standard DTQW evolution, identical coin oper- 
ation is used during each step of the walk making it a 
periodic evolution. This order can be broken by choos- 
ing a different coin operation for each step of the walk 
or by choosing a different coin operation at each position 
space. For simplicity, we will consider a different coin op- 
eration during each step. Evolution of DTQW using dis- 
ordered coin operation can be constructed by randomly 
choosing a quantum coin operator for each step from a 
set of operators in the U(2) group. That is, 
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FIG. 2: (Color online) Comparison of the probability distri- 
bution of the QW evolved by assigning quantum coin oper- 
ation in a U(2) group with random parameters in a differ- 
ent range to each step of the walk to that of the Hadamard 
walk, (a) Hadamard walk, which is the standard form of 
the DTQW with identical coin (0 = 7r/4), (b) the walk us- 
ing parameters assigned from £ {0,n/2} for B£ t e,(, 
(c) the walk using parameters assigned from £, £ £ {0,n/2} 
and 9i £ {0,7r/4} and (d) the walk using parameters as- 
signed from £, £ £ {0,7r/2} and 62 £ {tt/4, vr/2}. Localiza- 
tion is seen in the case of (d). The distribution is after 200 
steps of the walk on a particle with initial state at the origin 

l*i» s ) = 7f(|0)+i|l))(8|Vo). 



S(B ett e tXt ®l) ■ ■ ■ S(Bt: Mx ®t) ■ ■ ■ S(B 5o: e oXo ®l)\* m ) 

(6) 

with randomly chosen parameters £, 9, £ £ {0,7r/2} for 
each step. Though the coin parameters are randomly 
chosen for each step, the evolution is unitary and in- 
volves interference of amplitudes, and the effect is seen 
in probability distribution, see Fig. [TJ 

From the numerical evolution of 100 step QW shown in 
Fig. [U although interference effect is seen, we note that 
the variance of the distribution is much closer to the vari- 
ance of the CRW distribution. However, by restricting 



the range of the coin parameters that can be used for 
the walk, the probability distribution can be localized or 
made to diffuse in position space. One simple example we 
can consider is to by randomly pick different 9 for each 
step from a subset of a complete range of 9, with subsets 
61 £ {0,7r/4} and 9 2 £ {7r/4,7r/2}, for B Si g lX and B it e 2 x 
respectively, while other parameters are still picked from 
the complete range, £, £ £ {0,7r/2}. The probability dis- 
tribution obtained is shown in Fig. [2j The QW using 
9\ £ {0, 7r/4} diffuse in position space without any sharp 
peaks when compared to the standard QW evolution us- 
ing identical coin operation with two sharp peaks. The 
walk using O2 £ |7r/4, 7r/2|, however, localizes the dis- 
tribution around the origin [24| . For the numerical evo- 
lution, the parameter 9 for each random coin operation 
in the evolution was generated from a random number 
generator program with an equal probable appearance of 
any number in the specific range. 



IV. LOCALIZATION OF BOSE-EINSTEIN 
CONDENSATES 

Degenerate atomic gases have been used as a system to 
experimentally implement a number of basic models in 
condensed matter theory. The possibility to create both, 
ordered and disordered lattice potentials in higher phys- 
ical dimensions, the control of interatomic interactions, 
and the possibility to measure atomic density profile via 
direct ima ging are the key advantages of atomic quantum 
gases (cf. (Hi,[13)- 

Recently, atomic gases were used to test and implement 
various quantum information protocols; particularly, by 
using a single atom in an optical lattice, experimental im- 
plementation of DTQW has been reported [13] ■ Single 
laser cooled caesium (Cs) neutral atoms were determin- 
istically delocalized over the sites of a one-dimensional 
spin-dependent optical lattice. Initially, the atoms dis- 
tributed among the axial vibrational state was prepared 
in the |0) = \F = 4,ttif = 4) hyperfine state by opti- 
cal pumping, where F is the total angular momentum, 
and rriF is its projection onto the quantization axis along 
the dipole trap axis. The resonant microwave radiation, 
a 7r/2 pulse, was used to coherently couples this state 
to the |1) = \F = 3, nip = 3) state and initialize the 
system in the superposition (|0) + i\l))/y/2<S> \ipo)- The 
state-dependent shift operation is performed by contin- 
uous control of the trap polarization, moving the spin 
state |0) to the right and state |1) to the left adiabati- 
cally along the lattice axis. After t steps of coin operation 
and state-dependent shift, the final atom distribution is 
probed by fluorescence imaging. From these images, the 
exact lattice site of the atom after the walk is extracted 
and compared to the initial position of the atom. The 
final probability distribution to find an atom at position 
x after t steps is obtained from the distance each atom 
has walked by taking the ensemble average over several 
hundreds of identical realizations of the sequence. 
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FIG. 3: (Color online) Comparison of a standard QW evo- 
lution (Hadamard walk) and a localized QW (LQW) after 
100, 200, and 400 steps respectively. For a walk using ran- 
domly picked parameters 62 £ {7r/4, 7r/2} for each step of 
the walk, distribution remains localized near the origin irre- 
spective of the number of steps (time ). The distribution is 
after the walk on a particle with initial state at the origin 
\^ins) = ^75 (|0) ® IV'o)' Inset is a closer view of local- 

ized QW distribution between position -30 and +30 for 100, 
200 and 400 steps, respectively. 
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FIG. 4: (Color online) Localization length obtained by calcu- 
lating the ratio of the variance of the probability distribution 
LQW to that of the QW using an identical coin operation for 
each step. For each plot different 6 in the coin operation is 
picked for an identical coin operation. 



Experimental complexity aside, the preceding protocol 
can be adopted to other species of ultracold atoms and 
the non-interacting BEC as an initial state. For exam- 
ple, rubidium ( 87 Rb) atoms in an optical trap with state 
|F = l,m F = 1) = |0) and \F = 2,m F = 2) = |1) can 
be coherently coupled and moved in position space to 
implement DTQW. Implementing DTQW on the BEC 



involves evolving BEC into macroscopic superposition 
(Schrodinger cat) state [H-[6l| after each shift operation 
which delocalizes the state in position space [53| . Evolv- 
ing the BEC into a macroscopic superposition state to 
implement DTQW involves applying a rf pulse to the 
system to transfer the atoms part of the way between 
states |0) and |1). The duration of the pulse is kept 
much shorter than the self-dynamics of the condensate. 
This will only evolve each atom into the superposition of 
state |0) and |1), and the corresponding N atom quan- 
tum state is a product of single-particle superposition, 
that is, it is still a microscopic superposition. However, as 
this initial state evolves under the nonlinear Hamiltonian 
that governs the BEC with attractive interparticle inter- 
actions, it reaches a macroscopic superposition in which 
all atoms in a given lattice are simultaneously in level 
10) and |1>, \^ BEC ) = (l/V2)[|iV|o>,0|i)) + 10,0),^!))] 
|53l |6l| . Therefore, interatomic interactions are impor- 
tant to evolve the BEC into a macroscopic superposition 
state as discussed above. A DTQW on a noninteracting 
BEC as an initial state will implement the walk at an 
individual atom level and lose the coherence. We should 
note that the two experimental demonstration of local- 
ization [12, LH[ were done on a noninteracting BEC. 

In a scheme proposed to im ple ment a DTQW on ul- 
tracold atoms in a BEC state J53[, it is rst evolved to a 
macroscopic superposition state and a stimulated Raman 
kick, i.e., two selected levels of the atoms are coupled 
to the two modes of counterpropagating laser beams to 
coherently impart a translation of atoms in the position 
space. After each translation, the wave packet is again 
evolved into the macroscopic superposition state at each 
lattice x, where the number of atoms n x < N and the 
process is iterated to implement a large number of steps 
of QW. When n x is very small, the atoms will evolve to 
the superposition only at an individual atom level, i.e., 
microscopic superposition. With a certain modication 
to this scheme, that is, by evolving atoms to the super- 
position of the states at an individual atom level and 
implementing the shift operation before the interatomic 
interaction takes over to form a macroscopic superposi- 
tion, the QW at an individual atom level can be realized 

To observe localization, the periodic evolution using a 
7r/2 pulse as quantum coin operation during each step of 
the walk in the preceding protocols is broken by randomly 
picking the pulse from the range {7r/4, 7r/2}. Therefore, 
during each step, the coin operation evolves the state to 
an unequal superposition state such that the constructive 
interference effect and amplitude are directed toward the 
origin. In Fig. [3j we compare the QW evolution using an 
identical coin operation for each step (Hadamard walk) 
and localized QW (LQW) after 100, 200, and 400 steps 
respectively. For a walk using randomly picked parame- 
ters 62 € {7r/4,7r/2} for each step of the walk, distribu- 
tion remains localized near the origin irrespective of the 
number of steps ( time ) . 

In Fig. [5J the localization length L; oc is numerically 



■5 



obtained by taking the ratio of the variance of the prob- 
ability distribution LQW to that of the QW using iden- 
tical coin operation for each step, 



where 62 is the randomly picked quantum coin parameter 
from range {7r/4, it/ 2} for each step and 8 is the identical 
quantum coin parameter throughout the walk evolution. 



V. CONCLUSION 

We have discussed an approach to dynamic local- 
ization of ultracold atoms in a one-dimensional lattice 
under the inuence of DTQWusing disordered quantum 
coin operations. We introduced a DTQW model in 
which a random coin parameter is assigned to each step 
of the walk to break the periodicity during the walk 
evolution. By picking the coin operation from a different 
range of parameters, we have shown that the DTQW 
on a two-state particle in a one-dimensional lattice 
can be diffused or strongly localized in position space, 



respectively. We have shown that these behaviors of the 
DTQW can be efciently induced without introducing 
decoherence into the system. Using ultracold atoms in 
an optical lattice as a physical system, we have discussed 
implementation of a DTQW exibility in control over the 
experimental parameters to congure walk to diffuse or 
localize in position space. We have discussed implemen- 
tation of DTQW on atoms at an individual level in a 
BEC and on a BEC retaining the macroscopic coherence 
(BEC) state throughout the evolution. From this, we 
can conclude that the localization can be observed 
in a BEC with noninteracting atoms and interacting 
atoms, respectively. In general, the disordered coin 
operations (microwave pulses) in DTQW on atoms can 
be made to mimic the random media localizing the BEC. 
This approach can be adopted to any other physical 
system in which a controlled disordered DTQW can 
be implemented, broadening the spectrum of possible 
application of DTQW to study dynamics and phases in 
physical systems. 
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